PROCEEDINGS 

°F SCIENCE 



Pseudoscalar Flavor-Singlet Physics with Staggered 
Fermions 



UKQCD Collaboration 

Eric B. Gregory*! Alan Irving, Chris M. Richards 

Department of Mathematical Sciences, University of Liverpool, Liverpool, L69-7ZL, UK 
E-mail: |e . gregory@physics . gla . ac . uk 



Craig McNelle 

Department of Physics and Astronomy, University of Glasgow, Glasgow G12-8QQ, UK 

Allstair Hart 

School of Physics (JCMB), University of Edinburgh, Edinburgh. EH9-3JZ, UK 



Accurately calculating the mass of flavor-singlet meson states from numerical lattice simulations 
is an important milestone for lattice QCD. Careful measurement of the full pseudoscalar flavor- 
singlet propagator is also a crucial step in understanding the dynamics of the fermion sea on 
the lattice, in particular for potentially non-trivial formulations such as with 2h-1 -flavor staggered 
fermions. We briefly describe details of a dynamical QCD calculation using improved staggered 
fermions, with 30,000 trajectories, that was run for our studies of flavour singlet mesons. 



The XXV International Symposium on Lattice Field Theory 
July 30-4 August 2007 
Regensburg, Germany 



* Speaker. 

' University of Glasgow from October 2007 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Pseudoscalar Flavor-Singlets with Staggered Fermions 



Eric B. Gregory 



1. Introduction 

An first-principles calculation of the mass of the Tj' meson has long been a goal of lattice 
quantum-chromodynamics. Such a calculation would shed light on the dynamics of the fermionic 
sea and the mechanism by which fluctuations raise the mass of the rj' mass over that of the pion 

ii- 

A number of collaborations have performed simulations of pseudoscalar singlet systems with 



Nf = 2 flavors of dynamical fermions, e.g.[|3|, ^, |, |10[ [11|]. Recently the JLQCD/CP-PACS collab- 
oration reported on a preliminary calculation of the rj and rj' mesons with 2+1 flavours of Wilson 
fermions 1 12]. This article is based on the work with 2 + 1 flavors described more fully in [|l|]. 

There are compelling reasons to attempt such a calculation with improved 2+1 flavors of stag- 
gered fermions. Improved staggered fermions have an impressive track record of accurate hadron 



spectroscopy of flavour non-singlet quantities with light dynamical quarks [14, 15], suggesting it 
would be a good testing ground for the difficult task of singlet spectroscopy. The flavor-singlet 
propagator contains disconnected diagrams, which are inherently noisy in numerical simulations 
and require long Monte Carlo timeseries to measure with precision. Hence it is imperative that any 
program intent on singlet spectroscopy be capable of fast simulation of dynamical quarks. 

The staggered fermion program, however, is nagged by questions about the validity of the 
fourth-root trick which is employed to convert the native four tastes of staggered fermions to the 
desired Nf flavors. Recent theoretical work suggests that in the continuum limit [|l^, [l7|] problems 
with the fourth root may vanish, but if not, the pseudoscalar flavor-singlet system is a likely arena 



to find evidence of any inconsistencies [18, 1 



With this in mind we have endeavoured to simulate the tj/t]' system with Nf = 2 + 1 flavors 
of improved staggered fermions. In this article we report on the conclusions of the first stage of this 
project — an attempt at pseudoscalar flavor-singlet spectroscopy on the library of MILC "coarse" 
lattices, and discuss the prospects for the project in its second stage, namely measurement of the 
pseudoscalar flavor-singlet on a timeseries consisting of roughly 30000 trajectories of 2+ 1-flavor 
improved staggered fermion configurations. 

2. Theoretical background 

In general, the propagator for the pseudoscalar singlet meson with Nf flavors of quarks is given 

by 

Nf Nf 
Gsp{x',x) = {J^qi{x')r5qi{x') £ q j{x)r5qj{x)). (2.1) 

i=\ 7=1 

From the two types of contractions in this expression we get Nf connected diagrams: 



and Nj disconnected terms: 



(E ) ^5 qi ix)J^qj{x)Y5qjix)), (2.2) 
j 



(L qi {x! ) Ysqi {x ) £ qj (x) 75 qj (x)). (2.3) 
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In the Nf = 2 flavor-symmetric case, the connected term is identical with the pion propagator. 
So we can write: 

Gsp{x',x) = NfC{x',x) -NjD{x,x), (2.4) 
where the extra fermion loop in the disconnected diagram gives rise to the relative minus sign, and 

Gnp{x ,x) = NfC{x ,x). (2.5) 

A significant difference between the connected and disconnected terms is the role of the sea 
quarks bubbles in the disconnected term. One way to highlight the effects of the sea quarks is 
to take the ratio of the disconnected to connected terms. We expect that at large Euclidean time 
separations in full QCD 

GNp{t)=ANpe~"^'-' and Gsp{t) = Aspe'^'^' . (2.6) 

So we would expect the ratio of disconnected to connected parts to go as 

^t) = = 1 _ ^,-iMs.-M,.)t (2.7) 

NfC{t) Anp 

in the SU [Nf) flavor-symmetric case. In the relevant Nf = 2 + 1 case we would instead consider a 
modified ratio: 

By constructing this ratio we can highlight the behaviour of the sea quarks. For example, if 
instead of simulating full QCD we perform a quenched simulation, we should expect to find 

D(/,x) _ jml-aM^p) ml + aMlp 
^' C{x',x) 2Mmp 2M^p ' ^ ' 

where a is the parameter of the kinetic term of singlet pseudoscalar meson [|20|], and mo is the 
difference between the masses Mgp and M^p. 

By measuring the D/C ratio in (2^) we might expect to find any problems with the fermionic 
sea introduced by the fourth-root trick. 



3. Simulation details 



In the first phase of this project, we measured the pseudoscalar singlet connected correlator 
and the fermionic loops required to construct the pseudoscalar singlet disconnected correlator on 
Nf = 0,3 and 2+1 flavor lattices primarily from MILC's "coarse" ensemble library. The specific 
ensembles examined are listed in Table |l|. 

To more fully understand the fluctuations inherent in disconnected correlator measurements, 
we extended the quenched ensemble from 408 configurations to 6154 configurations. 

We examined the 75 1 taste-singlet state. This is generated with a four-link operator A^^^i 
which displaces the quark and anti-quark sources to opposite corners of the hypercube (and applies 
appropriate Kogut-Susskind phases). The parity partner of the 75 (8) 1 state is exotic so correlators do 
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10/^2 


L^xT 




arrival 


Nctg 







8.00 


20^ X 64 




0.020 


408 







8.00 


20^ X 64 




0.050 


408 - 


-^6154 


2 


7.20 


20^ X 64 


0.020 


0.020 


547 




2+1 


6.76 


20^ X 64 


0.007, 0.05 


0.007, 0.05 


422 




2+1 


6.76 


20^ X 64 


0.010, 0.05 


0.010, 0.05 


644 




2+1 


6.85 


20^ X 64 


0.05, 0.05 


0.05, 0.05 


369 





Table 1: Ensembles used for singlet calculations. 



not include contributions from an oscillating partner state. We use simple point sources to measure 
the connected correlators. 

To measure the disconnected correlator we use stochastic Gaussian volume sources. Following 
Venkataraman and Kilcup we use a noise-reduction trick that almost completely removes the 
component of the variance that is due to the stochastic sources. The Venkataraman Kilcup variance 
reduction (VKVR) trick involves replacing the pseudoscalar loop operator (T]Ayj0iM^^T]^),j with 
m{riAy^0iM^^M^^^ri^)rj. The staggered property of connecting only odd to even lattice sites 
ensures that the two have the same expectation value, but the latter has greatly decreased variance. 



4. Results 



An initial calculation of the D/C ratio (2^) on the 408 quenched j3 = 8.00 lattices with quark 
mass am = 0.05 was not convincingly linear (bursts in Fig. [T]). To explore why this was the case 
we extended the ensemble to 6154 configurations by generating more lattices in ten Monte Carlo 
streams and calculated the D/C ratio on the full ensemble (bold inverted triangles in Fig. [T]). With 



the increased statistics the D/C ratio is now consistent with the form on (2.9). 

We separately extracted the non-singlet 75 ® 1 ground state from the connected correlator and 
found Mnp = 0.5180(3). In doing a linear fit to the quenched D/C ratio, and rescaling by s/Nf = 
\/3, we find mo = 0.76GeV. 

A similar calculation using the R{t)su3 ratio form in (2^) with the 644 j3 = 6.76 configurations 
gives a{Msp —M^p) = 0.24(9). From the connected correlator we get aMpjp = 0.571(2). Using 
a = 0.125fm we extract M^p = 1280(142)MeV. 

For the j8 = 6.76, am = 0.01,0.05 ensemble we also calculate the effective masses diagonal 
and off-diagonal elements of the variational propagator 



G(Ar) 



C,,(AO - 2D,,(A0 -V2Bg,,{At) 



(4.1) 



-V2D,g{At) C,,(AO-D,,(AO 

These are shown in Fig. ^ along with effective masses of the light and strange connected corre- 
lators. The variational method relies on all elements of G{AT) having the same ground state. We 
see in Fig. ^. that the diagonal elements do share a common ground state, corresponding to about 
990MeV. The off-diagonal elements of G{AT) are mixed disconnected correlators, and the effec- 
tive mass corresponds to about 600MeV. We would expect the lowest stable state in this channel to 
be the rj at 548MeV. 
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We suspect that the explanation for this discrepancy hes in the hmited statistics available in 
these ensembles. The great difficulty in measuring R{t) and the singlet propagators is in the fluctu- 
ations in the disconnected correlator. Unlike many observables in lattice simulations, disconnected 
correlator, as a product of two roughly Gaussian-distributed quantities (the loop operators), has a 
distinctly non-Gaussian distribution, whose peak remains at zero. The mean value is determined by 
the asymmetry which, when the correlation is small, is almost entirely in the tails, many standard 
deviations away from the mean. An example is shown in Fig. ^. Our experience with the quenched 
ensembles shows that several hundred configurations is generally not sufficient to resolve the dis- 
connected correlators to do meaningful spectroscopy in configurations of this size. We did not find 
significant autocorrelations in any ensemble. 
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Figure 1: D/C ratio for 6154 j3 = 8.00 quenched lattices with valence quark mass am = 0.05, and for 11 
subsets of 400 configurations. We highlight five of these subsets, including the original MILC configurations, 
as disagreeing with the mean by more than one a. 

Figure 2: D/C ratio for coarse MILC ensembles. For the Nf = 3 and Nf = 2 + I ensembles Rsu3 is 
displayed. For the quenched configurations the single-flavour/? = D/C is plotted. 
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Figure 3: Effective masses for the (75 1) channel for the light and strange quarks for the /3 = 6.76, 
m=0.01/0.05 data set. 



Figure 4: Histogram of 393856 D{At = 2) measurements on 6154 quenched lattices. 
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arrival 
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2+1 


6.75 
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0.006, 0.030 


00.006, 0.030 


6000 



Table 2: Long run dynamical ensembles for singlet calculations. 



5. Long dynamical ensemble 

To attack the disconnected correlators forcefully, we have generated a A'y = 2 + 1 flavor en- 
semble of 5081 configurations at intervals of six RHMC trajectories. These 24^ x 64 configurations 
with j8 = 6.75 and am = 0.006,0.03 have been generated on the UKQCD collaboration's QCDOC 
machine. We chose the strange quark simulation mass of arris = 0.03 to coincide with values used 
by MILC in newer coarse ensembles 

We have completed preliminary spectroscopy on a portion of this ensemble and find aMj^ = 
0.17478(35) and aMp = 0.5298(17). A preliminary estimate of the static quark potential gives 
ro/£? = 3.810(51). 




Figure 5: Topological chai^ge time-series and histogram for about 17000 trajectories of the long dynamical 
run. 

We have measured the topological charge on the ensemble and the time-history shows good 
tunnelhng with an autocorrelation time in the range of 30 - 50 trajectories (Fig.|5|). 



6. Conclusions and Outlook 

We have completed the first phase of a project to explore singlet physics with staggered 
fermions. We have tested several methods of extracting singlet masses from connected and dis- 
connected correlator combinations. We again refer the interested reader to the extended analysis in 
[13], including analysis of the other ensembles in Table |l|. 

Comparing the original and the extended quenched ensemble highlights the difficulties of 
working with "small" ensembles, and the need for long dynamical-fermion ensembles to calcu- 
late disconnected diagrams for singlet physics. We hope to be able to resolve the tj/tj' spectrum in 
the ensemble described in Section ^ 
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